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Abstract 

This work discusses the approximations frequently used to calculate the contribution to pion 
production in proton-proton scattering near threshold from irreducible pion re-scattering. The 
reference result is obtained from the quantum mechanical reduction of the corresponding Feynman 
diagram to Time-Ordered Field Theory diagrams. We conclude that the Distorted Wave Born 
Approximation is quite adequate at threshold energies and even above. We discuss the choices for 
the energy of the exchanged pion, both for its propagator and for the ttN amplitude describing the 
re-scattering vertex. The calculation employs a physical model for nucleons and pions — pseudo- 
vector coupling for the ttNN vertex and realistic amplitudes for the ttN re-scattering and for the 
NN transitions in the initial and final states. 
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I. INTRODUCTION 



Recently, high quality data on meson production at threshold are available, and new 
measurements of relative small effects on those reactions, as charge symmetry breaking, are 
accomplished. This situation motivates us to revisit technical details and approximations 
currently done in the theoretical calculations of the matrix elements, independently of the 
use of chiral perturbation theory or meson theory with variable phenomenological content. 

Interest in meson production was spurred originally by the underprediction of the cross 
section of the pp — > ppn° reaction. In FigQthe diagrams for the impulse (FigJTJi)) and the 
irreducible pion re-scattering (Fig^)) contributions to that reaction are represented. In the 
pp — > ppir reaction, the isospin selection rule suppresses the dominant low-energy isovec- 
tor Weinberg- Tomosawa 7rN term of the irreducible pion re-scattering diagram. Addition- 
ally, a negative interference between the remaining impulse and isoscalar pion re-scattering 
terms Q, 0, may occur. Consequently, the empirical data for pp —> ppn° near threshold 
is only explained by phenomenological non-pionic, thus short-range, two-body mechanisms 
introduced in Ref.^y], or alternatively by an off-shell 

This last result, concerning the spectacular enhancement of the cross section by the off-shell 
behavior of the nN re-scattering, was not confirmed in Ref.Ql where a microscopic model 
for the nN re-scattering was considered. The high quality data for pp — > ppn° becomes 
therefore specially interesting since it allowed to establish the importance of the short-range 
mechanisms represented in FigCloQ, 0, 0, Q- We note that given the high momentum 
transfer at threshold, these short-range mechanisms, app earing also in the short-range com- 

. Il2|. play a decisive role. Also, for the 
same reason, the convergence of the chiral expansion is questioned or at least calls for a 
redefinition of the expansion variable, as explained in Ref. |l3l| . 

All mechanisms described above are derived from Feynman diagrams. Nevertheless, the 
evaluation of the corresponding matrix elements for the cross section proceeds through non- 
relativistic initial and final nucleonic wave functions, within the framework of the Distorted 
Wave Born Approximation (DWBA). Therefore the calculations apply a three-dimensional 
formulation in the loop integrals for the distortion, which is not obtained from the field- 
theoretical Feynman diagrams. Namely, the energy of the exchanged pion in the re-scattering 
operator (both in the tiN amplitude and in the exchanged pion propagator) has been treated 
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approximately and under different prescriptions in calculations performed till now. A theo- 
retical control of the energy for the loop integration embedding the non-relativistic reduction 
of the Feynman 7r-exchange diagram in the non-relativistic nucleon-nucleon wave functions 



is then still needed, as pointed out before in Refs. |14t Il5| . In this paper we deal with the 
isoscalar re-scattering term for the pp — > ppn° reaction, near threshold. Although in this 
case the re-scattering mechanism is indeed small, the amount of its interference with the 
impulse term depends quantitatively on the calculation method. We note, furthermore, that 
the pion isoscalar re-scattering term, which is energy dependent, increases away from thresh- 
old, and that for the charged pion production reactions the isovector term is important, and 
depends also on the exchanged pion energy. Thus the knowledge gained from the application 
discussed here to the pp —>■ ppn° reaction near threshold is useful for other applications. 
Specifically, this work investigates 

i) the validity of the traditionally employed DWBA approximation. We will use as ref- 
erence the result obtained from the decomposition of the Feynman diagram into Time- 
Ordered-Field-Theory (TOPT) diagrams, and realize how the last ones link naturally to an 
appropriate quantum-mechanical DWBA matrix element; 

ii) the choices for the exchanged pion energy which are unavoidable in the three- 
dimensional non-relativistic formalism underlying DWBA, when the link described in i) 
is not realized. They are, namely, the static approximation, the on-shell approximation in- 
.educed urst in RM.B (used in Ref.H and labeled later |l4j the U E — approximation) 
and the fixed kinematics approximation ]J extensively used. In these three approximations 
the pion energy is, respectively, taken as zero, the difference between the final and the ini- 
tial on-mass-shell energy of nucleon which emits it, and this difference calculated exactly at 
threshold; 

iii) the numerical importance of the three-body logarithmic singularities of the exact 
propagator of the exchanged pion, which are not present when the approximations mentioned 
in ii) are considered. 

This paper generalizes the work of Refs. jl^l Il^ which considers a toy model for scalar 
particles and interactions, does not include contributions from negative-energy states, treats 
the nucleons as distinguishable and therefore pion emission to proceed only from one nu- 
cleon. Within this toy-model, the DWBA amplitude for pion production threshold is clearly 
dominant over the other contributions and, for the final-state interaction case, only the fixed 
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kinematics approximation (for both propagator and vertex) leads to small deviations from 
the exact amplitude. Also, in the region between the thresholds for one-pion and two-pion 
production, the most used approximations for the pion propagator are found to produce the 



wrong energy dependence 



15] 



Our calculation employs a physical model for nucleons and pions and investigates how 
much of the features mentioned above survive in a more realistic calculation which uses 
a pseudo- vector coupling for the rcNN vertex; the %Pt itN amplitude and the Bonn B 
potential for the nucleon-nucleon interaction. Section |H] makes the decomposition of the 
distorted irreducible pion re-scattering Feynman diagram into the TOPT diagrams and con- 
nects the field-theoretical diagrams to the quantum-mechanical matrix elements of DWBA; 
Section ITTT1 presents the results and Section HVI a summary and conclusions. 

II. FROM THE FEYNMAN DIAGRAM TO DWBA 
A. Final-state interaction diagram 

1. The amplitude 

The Feynman diagram for the reaction pp — > ppn where the NN final-state interaction 
(FSI) proceeds through sigma exchange is represented in part a of FigEl After the nucleon 
negative-energy states are neglected, it corresponds to the amplitude 

M FSI = - f dAq ' V(C') - - m 

J (2vr) 4 Wo) Q' -E 2 -u 2 -ieE l -E 7T + Q> -Lu l + ie y) 
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Q' - E 2 + F 2 - u a + ie Q' - E 2 + F 2 + u a - ie Q' - u w + ie Q' + uj t - ie' 

where the exchanged pion has four-momentum Q' = ^Q' Q , q'\ The term V (Q'q) is a short- 
hand notation for the product of the -kNN vertex with the ttN re-scattering amplitude. The 
functions F\ i2 stand for the energy of the final nucleons, E\ i2 for the energy of the initial 
nucleons and a>i 2 , uj n , co a for the on- mass-shell energy values for the intermediate nucleons, 
exchanged pion and sigma mesons, respectively. In Appendix A details on these functions 
are given. All quantities are referred to the three-body center-of-mass frame of the nNN 
final state. 

In order to perform the energy integration in the variable Q' we implemented a partial 



fraction decomposition of the integrand, which isolates the poles of the pion propagator. 
Finally for the integration we closed the contour in the upper-half-plane. 

Because the partial fraction decomposition of the propagators in Eq.(JTJ) was done prior to 
the integration over the variable Q' Q , we have the following outcome which is independent of 
the choice of the contour of this integration: the only terms from the decomposition which 
do contribute to the integral correspond to the ones with only one pole, which happens to 
be the Q' — u % — a/ q" 2 + mf or the Q' = —uj n = — \fq 2 + pion poles. The other terms, 
with nucleon poles and/or sigma poles, together or not with pion poles, have all these poles 
located on the same half-plane and consequently their contribution vanish. 

We stress at this point that this method for the energy integration implies effectively 
that the tiN re-scattering amplitude is evaluated only for on-mass-shell pion energies. In 
this way, off-shell extrapolations which are not yet solidly constrained are avoided. Other 
methods may need the contribution of the off-shell amplitude in the integrand with the form 
shown in Eq.Q. But the net result is the same, provided that all the contributions from 
all (nucleon, sigma and pion) propagator poles are considered. So far, calculations did not 
consider the pion propagator poles, since they approximate that propagator by a form free 
of any singularity. Consequently, they exhibit a strong dependence on the nN amplitude at 
off-mass-shell energies of the incoming pion. 

After the Q' Q integration, one obtains for Eq.((TJ): 

\aFSI f d 3 q' 1 

M ™FT = -y 7-rs 4^ x ( 
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with E tot = 2E = Fi + F 2 + E w and E = E X = E 2 . 

This equation evidences that there are six contributions to the amplitude. These six 
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terms, originated by the four covariant propagators of the loop, can be interpreted as time- 
ordered diagrams. They are represented by diagrams a\ to a§ in FigJ21 This interpretation 
justifies the extra subscript label TOPT for the Mpsi amplitude in Eq.(J2J). 



2. Extraction of the effective pion propagator 

From the six terms in Eq.©, the first four terms (corresponding to diagrams a% to of 
Fig|2J) have the special feature that any cut through the intermediate state intersects only 
nucleon legs. Thus, they may be identified to the traditional DWBA amplitude for the final- 
state distortion. In contrast, in the last two diagrams a$ and of FigHJ any cut through 
the intermediate state cuts not only the nucleon legs, but also two exchanged particles in 



flight simultaneously. They are called the stretched boxes [14 1. 

Because of the identification of diagrams a± to with DWBA, we may collect the four 
first terms of Eq.(J2J) and obtain what we may call the exact expression for the DWBA 
amplitude: 

MDWBA = \S^f[ V M G "] {E^Ei-l^-^-u/™- (3) 

Here T stands for the transition-matrix of the final-state interaction. We have used both 
sigma exchange, (with m a = 550MeV) as in Ref.jl^ 

1 1 



yDWBA _ _J_ 



(4) 



(E tot — F 1 — E„ — u 2 — u) a ) (E tot — F 2 — E v — u)\— u a ) 
which makes Eq. (J3J) coincide exactly with Eq. (J2J , and also the T-matrix calculated from the 
Bonn B potential. In the derivation of the integrand of Eq.© the propagators for the two 
nucleons in the intermediate state originated fused into only one overall propagator in the 
non-relativistic form, 

1 

GNN = {Ek + Eb-K-vi-u*)' (5) 
The function V (toy) includes the contribution of the two pion poles u n and — toy corre- 
sponding to two different time-ordered diagrams, 

y _ V (toy) {E\ - E n - LOi - toy) + V (-^tt) (E 2 -U)2- ^tt) ^ 

toy 

where the function V (toy) is the product of the nN amplitude with the ttNN vertex. The 
kinematic factors ~ El ~ jE,r ~" 1 ~ a; ' r and jE2 ~^ 1 ~" 2 may be interpreted as form factors. 



For the realistic model used here^j we have 

g\ \ E^uJv EnU n -q' • 
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v K) = v (?, c) = - (a ■ i) 2c - (* - A) ^ _ C3 

P) 

where Qtt — <7tt) is the four-momentum of the emitted pion. 

In the derivation of Eq. (J3J) from Eq. (0) the function G n for the pion propagator turns to 
be exactly 

G * = + f ] 2 - [(^ - 5 - f ) - ^ - 2 ' (8) 

which gives the form of effective pion propagator appropriate for a DWBA final-state cal- 
culation and can be written as 

Gn = [u x + E„-E + w v ] [E-Ut-WnY (9) 
Thus, in what follows Eq.® is taken as the reference result, and allows to investigate the 

nn 

effect of the most common approximations hitherto used for the pion propagator 3, 6J. These 
approximations correspond to three different choices for the energy of the exchanged pion 
in the exchange process. In a non-relativistic framework particles are always on-mass-shell 
but intermediate states can be off-energy- shell. Considering, however, on-energy-shell nu- 
cleon states for the nucleons in the intermediate state of the distorted diagram, corresponds 
to assume that u± = 2E — E n — uj 2 . This assumption originates the so called on-shell 
approximation, introduced in Ref.jsJ where the propagator is given by 

= (E-ujlf-cu 2 (on — shell approximation). (10) 

This has the form of the Klein- Gordon propagator, with the pion energy given by E — u>2, 
i.e. the difference between the initial and final on-mass-shell energy of the nucleon emitting 
the pion. When furthermore one takes the kinematics as the one at threshold, E — lo\ — 
E — uj 2 = E — M = ?y = one originates the fixed kinematics approximation, where the 
propagator is taken as 

— ^ m h — 2 (fixed kinematics approximation). (11) 

Naturally, the deviation between the two approximations increases with energy. Finally, one 
may consider an instantaneous pion exchange process, by taking the energy of the pion to 
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be zero. This defines the static approximation, as 



G st = —-jr (static approximation). (12) 
It is part of the usual non-relativistic approximation for exchange Feynman diagrams. 



B. Initial-state interaction diagram 

1. The amplitude 

The corresponding Feynman diagram (part b of FigOJ), for the NN initial-state interac- 
tion (ISI) proceeding through sigma exchange, after the nucleon negative-energy states are 
neglected, generates the amplitude: 

A' 1 1 
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(13) 



Q' - E 2 + F 2 - uj a + is Q' Q - E 2 + F 2 + u a - is Q' - uj n + ie Q' + u n - is 

where we have used a notation analogous to one used for the final-state amplitude Eq.((T}. 

In order to perform the integration over the exchanged pion energy Q' Q , as before, a partial 
fraction decomposition to isolate the poles of the pion propagator was done. By closing the 
contour such that only the residues of the Q' = ±u T = ±a/?' 2 + vn\ poles contributes, one 
obtains the amplitude: 
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The six terms in Eq. (jl4j) are interpreted as contributions from time-ordered diagrams 
represented in Fig|31 b\ to be- This interpretation justifies the extra subscript label TOPT 
for the Ai ISI amplitude. Although Eq.Q and Eq. (jl4j) are formally alike, for the initial state 
an extra pole is present (besides that from the nucleons propagator), since it is energetically 
allowed for the exchanged pion to be on- mass-shell. All the singularities are decisive for 
the real and imaginary parts of Eqs.(j2J) and (|T4"|) and were handled numerically through 
subtraction methods (Appendix B). 



2. Extraction of the effective pion propagator 

Analogously to the final-state case, the first four terms in Eq. (jl4)) (diagrams b\ to 64 of 
FigOl where any cut of the intermediate state intersects only nucleon lines) are identified 
with the DWBA amplitude for the initial-state distortion, and the last two (diagrams 65 and 
b$ of FigOl with two exchanged particles in flight in any cut of the intermediate state) to 
the stretched boxes. In other words, the decomposition obtained in Eq. (jl4|) allows to write 
the exact or reference expression for the DWBA amplitude for the initial-state distortion, by 
collecting the four first terms, which have intermediate states without exchanged particles 
in flight, i.e., 

1 f d 3 q' 



MnwBA = -- J ^3 [ V K) G *J {Ei + E 2 -u 1 -u 2 ) Tn ^ (15) 
where 

f> , s V (u) n ) (E n + Ft - u) X - uj n ) + V (-uj n ) {F 2 -002- LJ n ) 

V{lo 7T ) = (16) 

and again V (u^) stands for the product of the tiN amplitude with the ttNN vertex (Eq.( |7j) ). 

In Eq. (ll5j) we made the replacement of the exchanged particle potential by the iViV 
transition-matrix and the two-nucleon propagator 

Gnn = (T? p r (17) 

(Ex + E 2 —U\- uj 2 ) 

is the non-relativistic global iViV propagator. As we did for the final-state interaction, 
we extracted from Eq. the effective pion propagator to be included in a DWBA-type 
calculation. It reads for the initial-state distortion, 

^ = [^ + ^ + f] 2 -[(f -^ + ^)-^] 2 ' (18) 
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This is the reference form for the pion propagator for the ISI case to be compared with the 
mostly used approximations in the literature, already discussed in detail for the FSI case: 

G^ k = hi — - (fixed kinematics approximation) 

G on = ^ _p 2 ^_ u 2 (on — shell approximation) (19) 
G st = —^2 (static approximation). 

We tested these approximations numerically, as we did for the corresponding ones in Eq.([12|) 
referring to the final state. In the next section we present the results obtained. 

III. RESULTS 

A. Stretched Boxes vs. DWBA 

In all calculations the NN and (NN)n channels considered refer to the transition 3 P ~~ * 
( 1 S'o) so, which is expected to be dominant. We found that the DWBA amplitude is clearly 
dominant over the stretched boxes in the realistic model considered, as the dotted line in 
FigEl documents for the FSI case. It is interesting to compare this result with the one of 
Ref. 14]. As shown in FiglH the stretched boxes amplitude is less than 1%, of the total 
amplitude and therefore is about 6 times more suppressed than in the dynamics of the toy 
model used in that reference. Replacing the nN amplitude from %Pt by a simple contact 
amplitude the stretched boxes amplitude makes those boxes slightly more important, still 
they do not exceed 4% of the DWBA amplitude (solid line in FigEJ. 

In terms of the cross section, the weight of the stretched boxes relatively to DWBA is 
even smaller, of the order of 0.02% at most. This is seen in FigEJwhere we compare, for three 
different values of the energy, the cross section obtained with only the DWBA contribution 
(&dwba), with the ones obtained with only the stretched boxes terms (cr stretched ). In both 
cases considered, V^f N + V*^ N and V^ N + contact, respectively left and right panel in 
Fig0 the cross section with the DWBA terms is clearly dominant, and in a more pronounced 
way when compared to the less realistic case of Ref.jlJ], not only at threshold but even for 
higher energies as 440 MeV. 

However, we note that the ratios and energy dependence of the amplitudes and of the 
cross sections are significantly influenced by the nN amplitude used in the calculation. The 
stretched boxes are seen to be amplified by the contact nN re-scattering amplitude, due to 
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an interplay between the nN and the NN amplitudes. Relatively to the more realistic 
amplitude, the contact tiN amplitude gives a larger weight to the low-momentum transfer. 
A realistic ttN amplitude satisfies chiral symmetry. This implies cutting small momenta and 
giving more weight to the region of large momentum transfer, which however in turn is cut 
by the nucleonic interactions. The difference between the two nN amplitudes is clearly seen 
by comparing the behavior of each curve on the left panel of FigEl with the corresponding 
curves on the right panel, for small values of the NN interaction cut-off. 

In the case of the initial-state amplitude, the stretched boxes amplitude is also much 
smaller than the DWBA amplitude for the two cases shown in FigJHl The cross sections 
obtained with only the stretched boxes terms were found to be less than 1.2% of the DWBA 
cross sections, even for laboratory energies up to T iah ~ 440MeV (see Fig[7|). 

The results presented justify the DWBA treatment for pion production. 

B. Approximations for the energy of the exchanged pion 

Considering the final NN state-interaction case, within the model with scalar particles 
and interactions solved in Refs. jl^l all the approximations for the pion energy taken at 
the nN re-scattering amplitude (absorbed in the function V of Eq.©) overestimate the cross 
section. We found that these conclusions still hold in the more realistic case of a pseudo- 
vector pion-nucleon coupling with the %Pt nN re-scattering amplitude. This is shown in 
Fig|H| We remark however that realistic nN and NN FSI distortion makes the deviation 
between exact and approximate results to be very small. Since the Bonn potential (as 
all OBE potentials) corresponds to a non-relativistic reduction where a specific method was 
applied to make it energy independent, we tested the effect on the results of an interaction in 
the final state corresponding to an energy-dependent sigma exchange. The results are on the 
panel a) of FigJHl The conclusions are about the same when we describe the nucleon-nucleon 
interaction by the (energy independent) Bonn B potential scattering T-matrix, instead of 
the sigma exchange. 

In particular the results for the reference calculation (heavy full line), the static approx- 
imation (dashed-dotted line), and the fixed kinematics prescription (dotted line) on FigjHl 
practically coincide. This means that the choice of energy for the pion propagator is not very 
decisive for the FSI amplitude, confirming that the usual low-energy static or instantaneous 
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approximation for exchange diagrams is justified in that case. 

We note however from FigJSjthat the cross section is dramatically enhanced if the "ad- 
hoc" fixed kinematics energy is considered at the same time for the re-scattering amplitude 
and the pion propagator. The results show indeed that the fixed kinematics (light full line 
in FigJSJ) and in particular the on-shell prescription (light dotted line with bullets), when 
taken both for the re-scattering vertex energy and pion propagator, are nor adequate above 
threshold, even at energies close to it, and they overestimate the cross section by a large 
factor. This finding is consistent with the results of Ref. jjj where the on-shell approximation 
for the energy at the ttN amplitude and the pion propagator was introduced for the first 
time. 

As for the initial-state distortion, our calculations checked the expected result that the 
cross section for the initial-state interaction is much smaller than the cross section for the 
final-state interaction. Still, regarding the approximations for the energy of the exchanged 
pion, the cross sections with the various pion propagator prescriptions defined in the previous 
section follow the trend observed for the final state amplitude. In particular, when the %Pt 
irN re-scattering amplitude is considered all the approximations for the propagator in the 
ISI amplitude overestimate the exact cross section, given by the heavy full line on FigO 
This happens for a nucleon-nucleon interaction given by the (energy independent) Bonn B 
potential (panel b) of FigEJ), or by the energy dependent sigma exchange (panel a) of Fig EI). 

When we compare the results for the initial state interaction with the results obtained for 
the final state interaction, we recognize however that the deviation between the approximate 
results and the exact one is much more significant for the ISI amplitude. The on-shell G on 
approximation (and also the fixed kinematics approximation) gives a larger cross section 
(panel b) in FigEJ). Its effect is much more pronounced than it was for the final-state 
distortion. This difference can be understood since the initial state interaction induces high 
off-shell energies in the intermediate nucleons. This enlarges the gap between the G on or G$ k 
and the G exact calculations, compared to what happens for the FSI case, where the nucleons 
emit the outcoming pion before their interaction happens. 

Also, what is behind the two approximations G on or G^ k being worse representations 
of the exact amplitude for the ISI case than for the FSI case, is the physics related to 
the logarithmic singularities of the exact pion propagator. These singularities are present 
only for the ISI amplitude, but are absent in the case of the approximate forms for that 
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propagator. Their effect is visualized in Fig^] which represents the imaginary part of that 
amplitude. On Fig^^we represent the ratio of the total cross sections (from the FSI, ISI 
and their interference amplitudes) for three approximations considered by the total exact 
cross section. The dominance of the FSI amplitude for the cross section makes the fixed- 
kinematics and the on-shell approximation very close at threshold, since for that amplitude 
the two approximations coincide, very near threshold. The conclusion is that the total cross 
section can be affected as much as 30% by the energy choices for the exchanged pion, and 
that the effect is largest close to threshold. 



IV. SUMMARY AND CONCLUSIONS 

We have generalized the work done in Refs. 3, Q on the study of the approximations of 
the quantum-mechanical calculation of irreducible pion re-scattering on the reaction pp — > 
pp-iv . We considered a physical model for nucleons and pions, combining a pseudo-vector 
coupling for the tiNN vertex, the %Pt nN re-scattering amplitude, and a realistic potential 
for the NN interaction in the final and initial state. We found that the effect of the usual 
choices for the pion energy in the nN re-scattering vertex, which is not fixed by a non- 
relativistic formalism, and thus has to be derived from the four-dimensional relativistic 
approach, can be significant. 

Both for the final- and initial-state interaction, our results show that 

i) the DWBA formalism is quite adequate at threshold and even at higher energies, since 
this part of the full amplitude is clearly seen to be dominant over the stretched boxes. 
This result confirms quantitatively the usual DWBA treatment for pion production. 
This is independent of the model for the nN re-scattering amplitude. Nevertheless, 
relatively to less realistic models, a chirally constrained amplitude reinforces even more 
the relative importance of the DWBA amplitude in the total amplitude; 

ii) the choice for the energy in the pion propagator is not very decisive in the amplitude 
for the final state distortion. However, the amplitude for the distortion of the initial 
state shows an important sensitivity to this choice. The approximate cross section can 
deviate as much as 30% from the exact result. 

iii) the choice for the energy of the exchanged pion in the re-scattering vertex, as pre- 
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scribed from an extraction from the Feynman diagram, is crucial for the cross section 
strength. The fixed kinematics and the on-shell prescriptions for the re-scattering 
vertex overestimate overwhelmingly the exact amplitude. 

We conclude that for evaluation of the irreducible pion re-scattering, one has to use 
the reference amplitudes given by Eqs.Q and (JT5j) . The latter were obtained as non- 
relativistic reductions of the corresponding Feynman diagrams for the final- and initial- 
state interactions, respectively. In the future, the reference amplitudes corresponding to the 
resonance contributions for the charged n or the r\ meson production reactions, should be 
analyzed in the same way. 
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APPENDIX A: KINEMATICS 



For the final state NN interaction, the expressions for the on-shell energies of the particles 



are: 



L0 n 



ml + 



-p + f + q k 



M 2 + 



Qk =F 



ml + \q k - q u \ 



exchanged 7r 

intermediate nucleons 

a exchanged 
final nucleons 



F 1)2 = M + ^ 

For the initial NN interaction, the corresponding expressions are 



= 

^1,2 

UJrr = 



ml + 



q u + 



qk 



M 2 + |^| 



m l + I V - Qk\ 



F h2 = M+^ 



exchanged tt 

intermediate nucleons 
a exchanged 
final nucleons 



(Al) 



(A2) 



p and q^ are, respectively, the initial nucleon tri-momentum and the emitted pion tri- 
momentum. (q u ) is the relative tri-momentum of the two intermediate (final) nucleons 
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and E n = y m\ + \ q^f is the energy of the emitted pion. All quantities are referred to the 
three-body center-of-mass frame of the nNN final state. 



APPENDIX B: TREATMENT OF THE LOGARITHMIC SINGULARITY OF 
THE PION PROPAGATOR (ISI) 

A partial wave decomposition involving two angles arccos^i = Z (q u — qk,q^) an d 
arccosx 2 = Z(q^,qk) was implemented. In terms of X\ and x 2 the pion propagator for 
the FSI amplitude given by Eq.© reads 

g,=-( VA 1 1 ) < B1 > 

\qirh [X2) yi-XlJ \qirh {X2) V2-XiJ 

where a,\ = E — Fl ~ F2 — u 2 — a 2 = E + Fl ~ F2 — ojx + (ai and a 2 must be positive 
or zero) and I 2 (X2) = \q~^ — ~q~k \ = \/ q\ — ^QuQkX2 + q\- The roots of the denominator of 



are given by: 



2 rl2 2 

Vi = a \~'\ with P 2 = ml + % + I 2 {x 2 f and i = 1, 2 (B2) 

qnh [X2j 4 



Once G n is in this form, we have to deal with its poles at x\ = yi- We used the subtraction 
technique exposing the three-body logarithmic singularity: 

f fjB^A P[(xi)dXl = PV f l fM-f(y^) P . w+2 /(K , I2)Q Lfa) _ 

J-i yi - xi J -1 yi - xi 

-mP L { yx )f{y h x 2 ) (B3) 

where and are the Legendre polynomials of order L and Ql are the Legendre functions 
of the second kind of order L. These last functions exhibit logarithmic singularities which 
are given by the condition yi = ±1. We determine its solution(s) analytically. Then, the 
integral over the moving logarithmic singularities is handled by a variable mesh. When the 
number of singularities (always between zero and two) is n s , the interval of integration is 
divided into 2n s + 1 regions with breaking points given by the found singularities. In each 
one of the regions one considers a Gaussian mesh. 
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FIGURE CAPTIONS 



Fig 1 - Impulse, re-scattering and short-range processes which contribute to n 
production. 

Fig 2 — Decomposition of the Feynman diagram in terms of six time-ordered diagrams 
for the final-state interaction. The pion(sigma) field is represented by a dashed(solid 
double) line. The nucleons are represented by solid lines. The DWBA amplitude may be 
identified to the first four time-ordered diagrams (ai to a 4 ). The last two diagrams (a 5 to 
ae) are usually called stretched boxes. 

Fig 3 — Decomposition of the Feynman diagram in terms of six time-ordered diagrams 
for the initial-state interaction. The pion(sigma) field is represented by a dashed(solid dou- 
ble) line. The nucleons are represented by solid lines. The DWBA amplitude corresponds 
to the first four time-ordered diagrams (bi to 6 4 ) and the stretched boxes to the last two (b 5 
to b e ). 

Fig 4 — Importance of the stretched boxes compared to the DWBA FSI amplitude 
M.^c^pti as a function of the mass of the scalar particle for the final NN interaction. The 
ttN amplitude is contact re-scattering vertex (solid line) and the %Pt (doted line). The 
energy is taken at the pion production threshold. Absolute values of the amplitudes are 
considered. 

Fig 5 — Importance of the stretched boxes compared to the DWBA (FSI) for the cross 
section as a function of the mass of the scalar particle for the final NN interaction. The 
ttN amplitude is the %Pt amplitude (left) and a contact re-scattering vertex (right). 

Fig 6 — The same as FigHbut for the DWBA ISI amplitude M^opt- 

Fig 7 — The same as FigEl but for the ISI case. 
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Fig 8 — Effect of the choices for the energy of the exchanged pion when the %Pt rcN 
re-scattering amplitude is considered, as a function of the laboratory kinetic energy. The 
nucleon-nucleon interaction is sigma exchange (panel a)) and the Bonn B potential (panel 
b)). The cross section curves shown correspond to FSI alone. The heavy full line is the ref- 
erence calculation (Eq.®). The dashed-dotted line is the static approximation. The dotted 
line is the fixed kinematics prescription and the dashed line is the on-shell prescription. 
The light full line and the light dotted line with bullets correspond respectively to the last 
two prescriptions, taken not only for the propagator but also for the irN re-scattering vertex. 

Fig 9 — The same as FigJSJbut for the ISI amplitude. 

Fig. 10 — Imaginary part of the ISI amplitude AAdwba as a function of the dimen- 
sionless pion momentum. The exact result (full line) and fixed kinematics (dashed line), 
on-shell (short-dashed line) and static (dotted line) approximations are shown. 

Fig 11 — Ratio between total approximated and exact cross section as a function of the 
dimensionless pion momentum. The fixed kinematics (full line), on-shell (dashed line )and 
static (dotted line) cases are shown. 
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FIGURES 
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a) b) c) 



FIG. 1: Impulse, re-scattering and short-range processes which contribute to n production. 
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°1 0,2 03 04 05 a§ 



FIG. 2: Decomposition of the Feynman diagram in terms of six time-ordered diagrams for the 
final-state interaction. The pion(sigma) field is represented by a dashed(solid double) line. The 
nucleons are represented by solid lines. The DWBA amplitude may be identified to the first four 
time-ordered diagrams {a\ to 04). The last two diagrams (05 to a§) are usually called stretched 
boxes. 
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b 



b l h & 3 64 & 5 



FIG. 3: Decomposition of the Feynman diagram in terms of six time-ordered diagrams for the 
initial-state interaction. The pion(sigma) field is represented by a dashed(solid double) line. The 
nucleons are represented by solid lines. The DWBA amplitude corresponds to the first four time- 
ordered diagrams (b\ to 64) and the stretched boxes to the last two (65 to be). 
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FIG. 4: Importance of the stretched boxes compared to the DWBA FSI amplitude MfpQ PT , as 
a function of the mass of the scalar particle for the final NN interaction. The ttN amplitude is 
contact re-scattering vertex (solid line) and the %Pt (doted line). The energy is taken at the pion 
production threshold. Absolute values of the amplitudes are considered. 
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FIG. 5: Importance of the stretched boxes compared to the DWBA (FSI) for the cross section as 
a function of the mass of the scalar particle for the final NN interaction. The ttN amplitude is 
the %Pt amplitude (left) and a contact re-scattering vertex (right). 



24 




25 




FIG. 7: The same as Fig|S]but for the ISI case. 
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FIG. 8: Effect of the choices for the energy of the exchanged pion when the %Pt irN re-scattering 
amplitude is considered, as a function of the laboratory kinetic energy. The nucleon-nucleon 
interaction is sigma exchange (panel a)) and the Bonn B potential (panel b)). The cross section 
curves shown correspond to FSI alone. The heavy full line is the reference calculation (Eq.(j2J)). The 
dashed-dotted line is the static approximation. The dotted line is the fixed kinematics prescription 
and the dashed line is the on-shell prescription. The light full line and the light dotted line with 
bullets correspond respectively to the last two prescriptions, taken not only for the propagator but 
also for the irN re-scattering vertex. 
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FIG. 9: The same as Fig|S]but for the ISI amplitude. 
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FIG. 10: Imaginary part of the ISI amplitude Mdwba as a function of the dimensionless pion 
momentum. The exact result (full line) and fixed kinematics (dashed line), on-shell (short-dashed 
line) and static (dotted line) approximations are shown. 
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FIG. 11: Ratio between total approximated and exact cross section as a function of the dimen- 
sionless pion momentum. The fixed kinematics (full line), on-shell (dashed line) and static (dotted 
line) cases are shown. 
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